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ABSTRACT 

This paper gives an overview of the attempts to determine the distribution of dark matter in low surface 
brightness disk and gas-rich dwarf galaxies, both through observations and computer simulations. Observa- 
tions seem to indicate an approximately constant dark matter density in the inner parts of galaxies, while 
cosmological computer simulations indicate a steep power-law-like behaviour. This difference has become 
known as the "core/cusp problem", and remains one of the unsolved problems in small-scale cosmology. 
Subject headings: 



1. INTRODUCTION 

Dark matter is one of the main ingredients of the universe. 
Early optical measurements of the rotation of spiral galaxies 
indicated the possible pres ence of large amo unts of dark mat- 
ter in their outer parts (e.g. Ru bin et alJl9 78). though in many 
cases the r otation curve could also be explained by the stars 
alone (e.g. lKalnajslll983b lKentlll986h . Observations at even 
larger distances from the galaxy centers, using the 21 -cm line 
of neu tral hydrogen , definit ively confirmed the mass discrep- 
ancy dBosmal 1 19781 Il981alfbh . For an extensive review see 
ISofue & Rubin (2001). Most of these early observations con- 
centrated on late-type disk galaxies, which all share the prop- 
erty of having an almost constant rotation velocity in their 
outer parts (the so-called "fiat rotation curve"). As the dy- 
namical contribution of the stars and the gas is insufficient to 
explain the high rotation velocities in the outer parts, this im- 
plies that most of the observed rotation there must be due to 
some other material, the "dark matter". The observed con- 
stant velocity suggests that the dark matter in the outer parts 
of galaxies has a mass density profile closely resembling that 
of an isothermal sphere, i.e., p ~ r~ 2 . 

In the inner parts of galaxies stars are obviously present and 
they must be the cause of a (possibly large) fraction of the 
observed rotation velocity. This therefore leads to a transition 
from the inner parts where the stars contribute to (and in many 
cases dominate) the dynamic s, to the outer parts where the 
dark matter is important (e.g .,[Kalnais 1 983llKentll986lll987t 
Ivan Albada & Sancisi|[l986D . 

The rotation velocity associated with dark matter in the 
inner parts of disk galaxies is found to rise approximately 
linearly with radius. This solid-body behaviour can be in- 
terpreted as indicating the presence of a central core in the 
dark matter distribution, spanning a significant fraction of 
the optical disk. Some authors adopt a non-singular isother- 
mal sphere to describe this kind of da rk matter mass dis- 
tribution (e.g., Athanass oulaetal] 119871). while others pre- 
fer a pseudo- isothermal sphere (e.g., |Begeman et alj 1199 It 
Broeil sj | 1992|). Both mode ls describe the data well (see 
Kormendy & Freeman 2004), and by the late 1980s they had 
become the de facto description of the distribution of dark 
matter in (gas-rich, late-type) dwarfs and disk galaxies. 

In this paper we use the pseudo-isothermal (PI) sphere to 
represent the cored models, though this particular choice does 
not affect any of the discussion in this paper. The mass density 



distribution of the PI sphere is given by: 



PPi(r) = 



l + (r/R c )2' 



(1) 



where p (l is the central density, and Rq is the core radius of the 
halo. This density distribution leads to an asymptotic flat ve- 
locity Vqo given by = (A-nGpoR^)^ 2 , where G is the grav- 
itational constant. 

In the early 1990s, the first results of numerical iV-body 
simulations of dark matter halos based on the collisionless 
cold dark matter (CDM) prescription became available. These 
did not show the observed core-like behaviour in their in- 
ner parts, but were better described by a steep power-law 
mass-density distribution, the so-called cusp. Fits to the 
mass-distributions as derived from these early simu lations 
(IDubinski & Carlberdl99UINavarro et alJl996blfT99l indi- 
cated an inner distribution p ~ r a with a = — 1. (In the follow- 
ing we will use a to indicate the inner mass density power-law 
slope.) 

The results from these and later simulations are based on 
the (A)CDM paradigm, where most of the mass-energy of our 
universe consists of collisionless CDM in combination with a 
cosmological constant A. This ACDM paradigm provides a 
comprehensive description of the universe at large scales (as 
shown most recently by the Wilkinson Microwave Anisotropy 
Probe (WMAP) results; see ISpergel et all 120071) . However, 
despite these great successes, it should be kept in mind that the 
cusp and the central dark matter distribution are not predicted 
from first principles by ACDM. Rather, these properties are 
derived from analytical fits made to dark-matter-only numer- 
ical simulations. While the quality and quantity of these sim- 
ulations has improved by orders of magnitude over the years, 
there is as yet no "cosmological theory" that explains and cor- 
rectly predicts the distribution of dark matter in galaxies from 
first principles. 

The value a = -1 found in the early CDM simulations is 
very different from that expected in the PI model, where the 
constant density core (p ~ r ) implies a = 0. These two cases 
thus lead to two very different descriptions of the dark matter 
distribution in galaxies. The "cusp" (a = -1) models gives rise 
to a rapidly increasing "spiky" dark matter density towards the 
center, while the "core" model (a = 0) has an approximately 
constant dark matter density. The cusp model therefore has a 
rotation curve that will rise as the square-root of the radius, 
while the core model rotation curve rises in a linear fashion. 
The difference in shapes between the rotation curves of both 
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models is quite pronounced, and, in principle, it should there- 
fore be possible to identify CDM haloes in real galaxies by 
measuring their rotation curves. 

Over the last 15 years or so, much effort has been put into 
determining the central mass distribution in galaxies using 
their rotation curves, and comparing them with the outcomes 
of ever more sophisticated numerical simulations. To first or- 
der, one can summarize this work as observational determina- 
tions yielding slopes a ~ 0, while simulations produce a ~— 1 
slopes. This persistent difference is known as the "core/cusp 
controversy", sometimes also described as "the small-scale 
crisis in cosmology". The attempts to reconcile the observa- 
tions and simulations, either by trying to improve them, or 
by trying to quantify systematic effects or missing physics, 
are the subjects of this paper. I give a brief overview of past 
and present work dealing with the determination of the central 
dark matter density distribution in galaxies, with an emphasis 
on the observational efforts. An overview like this, touch- 
ing on many different topics in galaxy evolution, cosmology 
and computational astrophysics, is never complete, and only a 
small (but hopefully somewhat representative) fraction of the 
many papers relevant to this topic can be referred to in the 
limited space available. The rest of this paper is organised as 
follows: Sect. 2 gives a description of the results that numeri- 
cal simulations have produced over the years. Section 3 deals 
with the observational determinations of the dark matter den- 
sity distribution. Section 4 discusses physical scenarios that 
have been proposed to reconcile the core and cusp distribu- 
tions. Section 5 briefly summarizes the work discussed. 

2. COLD DARK MATTER CUSPS 

The presence of a cusp in the centers of CDM halos is one 
of the earliest and strong est results derived from cosm olog- 
ical N-body simulations. iDubinski & Carlbergl (11991b were 
among the first to investigate the density profiles of CDM ha- 
los and found that the inner parts of these simulated halos 
could be characterized by a power-law slope a = — 1. They 
did not rule out the existence of central cores, but noted that 
these would have to be smaller than the resolution of their 
simulations (~ 1.4 kpc). Subsequent simulations, at higher 
and higher resolutions, made the presence of cores in simu- 
lated CDM haloes incre asingly unlikely . 

A systematic study by Nava rro et al.l d!996bl 11997b of sim- 
ulated CDM halos, derived assuming many different sets of 
cosmological parameters, found that the innermost dark mat- 
ter mass density distribution could be well described by a 
characteristic a = -1 inner slope for all simulated halos, inde- 
pendent of mass, size or cosmology. A similar general result 
was fo und for the outer mas s profile, with a steeper slope of 
a = -3. iNavarro et al.l (119971) called this the "universal density 
profile" and it is described by 
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where p, is related to the density of the universe at the time 
of the time of halo collapse and R s is the characteristic radius 
of the halo. This kind of profile is also known as the "NFW 
profile". 

The corresponding rotation curve is given by 



V(r) = 



' ln( 1 + cx) — cx/(l + cx) 



(3) 



xpn(l+c)-c/(l+c)] ' 
with x = r//?2oo- This curve curve is parameterized by a ra 



dius Z?200 an d a concentration parameter c = Rioo/Rs- Here 
7?2oo is the radius at which the density contrast with respect 
to the critical density of the universe exceeds 2 00, roughly the 
virial r adius; V200 is the circular velocity at R200 (INavarro et al.1 
1996b). The parameters c and V200 are tightly related through 
the assumed cosmology. Indeed, one can be expressed as a 
funct ion of the other, with only a small scatter dBullock et al.l 
12001 . That is, the range of (c,V2oo) combinations that de- 
scribes "real" CDM rotation curves is tightly constrained by 

the ACDM cosmology. 

Simulations by [Moore et al. (1999) indicated an even 
steeper inner slope. They found that their simulated halos 
could be best described by a function 



pym(r) = 



Pi 



(r/R s YHl + r/R s y- 



(4) 



i.e., with an inner slope a = -1 .5 and an outer slope a = -3. 

The difference between these two results indicated that is- 
sues such as numerical convergence, initial conditions, anal- 
ysis or interpretation could still play a role in defining the in- 
ner slope. As ever more powerful computers and increasingly 
higher resolution simulations became available, the value 
and behavior of the inner slope of CDM halos has there- 
fore been extensively discussed in the literature. For exam- 
ple, to give but an incomplete listing of the many papers 
that have appeared on this topic, Klvpin et al] d2001h derived 
slopes a = -1.5 for their simulated ha l os. Fr om phase-space 
density arguments, Tayl or & Navarro! ([2001) argue that the 
density profile should resemble an NFW profile, but con- 
vergin g to an inne r slope a = -0.75, instead of the a = -1 
value. iCohn etafl ([2004) investigated low-mass haloes and 
found that they were bes t desc ribed using NFW profiles (i.e., 
a = -1). iDiemand etail d2005l) found that CDM halos have 
cusps with a slope a ^ — 1.2. 

Many studies assumed that the central cusp consisted of 
a region where the mass density behave d as a power-law 
with a constant slope. INavarro et al.l (120041) and Hayas hiet alJ 
(2004) suggested that this did not have to be the case. They 
did not find evidence for an asymptotic power-law slope, 
but instead noted that the slope kept getting shallower to- 
wards smaller radii without converging to a single asymptotic 
value. At the smallest resolved radii they derive slopes of 
~ —1.2 for "galaxy-sized" halos (as measured at ~ 1.3 kpc), 
and ~ -1.35 for "dwarf galaxy" halos (as measured at ~ 0.4 
kpc). These values are significantl y steeper than the o riginal 
NFW slope, but not as steep as the Moore et al. ( 1999) value. 
INavarro et alJ ([2004) introduce a new fitting formula to quan- 
tify their results. For reasonable choices of its input parame- 
ters, this formula yields an extrapolated slope of a ~ —0.7 at 
r ~ 0.01 kpc. 

Stoehr (2006) also finds a gradual turn-over in slope to- 
wards smaller radii. Though his simulations formally resolve 
only radii ~ 1 kpc (where a slope of a ~ — 1 is measured), an 
extrapolation of his favoured fitting function towards smaller 
radii results in a decreasing slope ending up as a flat slope 
( a = 0) around r ~ 0.0 1 kpc. 

iMerritt et ail d2005l) and iGraham et ail d2006l) showed tha t 
the density distribution presented in INavarro et alJ (2004) 
could be equally well described by a Sersic function. In the 
context of CDM halos they refer to this function as an Einasto 
model. For completeness, this profile is given by 



PEin(r) = p e exp (-d„ (r/r e ) 1/n - 1 



(5) 
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where n determines the shape of the profile, and d„ is a func- 
tion of n which enables the use of the density p e measured 
at the effective radius r e . The latter is defined as the radius 
of the volume containing half of the total mass. In terms of 
observationally more accessible quantities this can be written 
as 

pm*(r) = P-2,Ein exp {-In (r/ r- 2tEm ) 1/ " ~ 1 ) , (6) 

where r_2,Ein is the radius at which the logarithmic derivative 
of the density profile equals -2, and p_2,Ein is the density at 
that radius. The two versions of radius and density are related 
by yO-2,Ein = p e exp(d„-2n) and r_ 2 . E in = (2n/d„)"r e . 

A further discussion of this model is beyond the scope of 
this paper, except to note that for typical parameterisations of 
CDM galaxy halos one derives a slope of a <~ —1.3 ± 0.2 at a 
radius of 1 kpc and of a ~ -0.9 ± 0.2 at a radius of 0. 1 kpc. 
The precise values depend on the exact values of n and r-i\ 
the value s just listed assume 4 < n < 8 and r_2 = 10 kpc, as 
shown in Graham et al. (200^) 1 . 

Amongst the highest resolution measureme nts of the inner 
slop e of CDM halos so far are those by Nav arro et al.l (2008) 
and IStadel et all (12008b. The former were d one as part of 
the Ac quarius project (jSpringel et al. 2008). iNavarro et al.l 
(2008) found no convergence to a single asymptotic slope. 
Rather, as before, the slope keeps decreasing with decreasing 
radius. At r = 0.1 kpc, which is approximately the smallest 
reliably resolved radius, the slope has a value a ~ 0.85. At 
r = 1 kpc the value is a ~ 1 .4. An Einasto profile with n = 5.9 
provides a good fit to the change in slope with radius. 

In an indepen dent, but equally detailed simulation, 
IStadel et"afl (120081) find a similar behavior, as well as com- 
parable slope values. They quote a slope a = -0.8 at 120 pc, 
and a = -1 .4 at 2 kpc. 

Even though the details of the simulations, the analytical 
fits and the interpretation and analysis differ, we can still draw 
some conclusions from the above discussion. All simulations 
and fitting functions considered here produce slopes a < — 1 at 
a radius of 1 kpc. At radii less than 1 kpc the most recent sim- 
ulations tend to produce slightly more shallow slopes where a 
typical value seems to be a ~ -0.8 at 0.1 kpc. From an ob- 
server's perspective, all models described here produce slopes 
a ~ -0.8 or steeper at the smallest observationally accessible 
radii, and will thus all produce results very similar to those 
derived using a "standard" NFW profile. In the simulations, 
radii less than 0.1 kpc cannot yet be reliably resolved, and the 
values of the slope derived there depend on the validity of the 
assumed analytical fitting function. 

3. OBSERVATIONS 

3.1. Early measurements 

The first comparisons of the HI rotation curves of gas-rich 
dwarf galax ies with those pre dicted by CDM profiles were 
presented in lMoord(fl994l) and lFlores & Primacld(ll994l) . The 
dynamics of these galaxies are dominated by dark matter, and 
they are therefore thought to be good probes of its distribution. 
Both studies note a large discrepancy between the observed 
rotation velocities and those predicted, especially in the inner 
parts. They show that the PI model gives a superior descrip- 
tion, implying that the halos of these late-type dwarf galaxies 
are best characterized by an approximately constant-density 

1 Less steep slope values listed in the same paper are derived assuming 
r_2 = 100 kpc and are thus more appropriate for giant or group-sized CDM 
halos 



core. iMoord ([1994) briefly addresses some of the observa- 
tional uncertainties that might affect the data, such as reso- 
lution, projection effects due to inclination, and the effects 
of pressure support, and concludes they are not significant 
enough to affect the results. He also notes that it is conceiv- 
able that, during the galaxy formation process, gas settling in 
the halo will have affected the dark matter distribution. Usu- 
ally this is thought to take place in the form of a pr ocess called 
"adiabatic contraction" (Blu menthal et al.| [l986). which has 
the effect of contracting the inner dark matter distribution (in- 
creasing the density). If currently derived halo properties are 
the result of this process, then the initial halos must have been 
o f even lower density, ex acerbating the discrepancy. 

INavarro et al.l (Il996al) argue that baryonic processes might 
be the cause of the observed core distribution. They use N- 
body simulations to model the effect of star formation on the 
baryons and the dark matter, and find that a central dark mat- 
ter core can be created if a large fraction of the baryons is 
suddenly expelled into the halo. They estimate that star for- 
mation rates of up to 10 M Q yr -1 are needed over a dynamical 
timescale of a galaxy for the process to have the desired im- 
pact. 

Afte r analysing t he same four d warf galaxies that Moore 
(1994) investigated, Burkert (1995) comes to the conclusion 
that their rotatio n curves, after appropriat e scaling, are self- 
similar (see also Salucci & Burkert 2000). He notes that it 
is unlikely that baryonic blow-outs and mass-flows can cause 
this kind of behaviour, unless fine-tuned, and attributes the 
slow rise of the rotation curve to the intrinsic properties of 
the dark matter (i.e., a dark matter core). He also finds tenta- 
tive evidence that the mass density in the outer parts of these 
dwarfs drops off as r" 3 (consistent with CDM), and not with 
r~ 2 (as suggested by the asymptotically flat rotation curves of 
spiral galaxies). Based on this, he introduces what has since 
become known as the "Burkert-profile": 



PBui- - 
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(7) 



with pa the central density, and ro the scale radius, similar 
to the core radius R c of the PI model. This model is thus 
characterised by a = in the inner parts, and a = -3 in the 
outer parts. 

A paper by Gelato & Sommer-Larser] (119991) presents a 
more detailed analysis of the baryonic blow-out process and 
its impact on the halo structure. They attempt to reproduce the 
observed rotation curve of DDO 154 by simulating NFW ha- 
los, and subjecting them to the effect of violent gas outflows, 
which are simulated by suddenly changing the disk potential 
in the simulations. In order for the final rotation curve to re- 
semble the observed curve, they need to suddenly blow out 
33 to 75 percent of the initial disk material. They note that 
they can reproduce the rotation curve of DDO 154 for a rather 
wide range of blow-out scenario s, and suggest that the fine- 
tuning argument put forward by Burkerjl (119951) may not be 
applicable. 

The blow-out process immediately gives rise to a number 
of observational consequences. Firstly, a period of star for- 
mation intense enough to blow out the majority of the baryons 
should leave behind a substantial (by now) old stellar popula- 
tion. In practice, however, these dwarfs are seen to be dom- 
inated by a young stellar population. Is there a discrepancy 
here? Secondly, what happens to the baryons that are blown 
out? Do they stay in the halo? Presumably they will be in the 
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form of hot gas. Is this hot gas detectable? Or, if the gas cools 
down, can we see it raining back on the disk? Discussing 
these questions in detail is beyond the scope of this paper, 
but note that there are a number of starbursting dwarf galax- 
ies in our local universe where one can attempt to study these 
phenomena directly. As an example, lOtt et al.l ([2005) anal- 
yse the properties of the hot and cold gas in 8 dwarf galaxies 
that are in a starburst phase. They show that outflows of hot 
gas are possible, but also that the presence of (tidal) cold gas 
can again confine the hot gas. In the galaxies in their sam- 
ple, approximately 1 percent of the total ISM is in the form 
of hot, coronal gas. The outflows are found to be efficient in 
removing hot, metal-rich gas. Whether these processes can, 
in the early universe, also remove the bulk of the ISM is still 
an open question. It is clear that with this kind of analysis we 
are no longer in the realm of cosmology, but are dealing with 
"messy" astrophysics. 

Fortunately, there is an alternat i ve way to study the dark 
matter distribution. Navarro et al. (1996a) already note that 
the blow-out process can only be effective in dwarf galax- 
ies. In more massive galaxies, such as spiral galaxies, the 
potential well is too deep to efficiently remove the gas. Find- 
ing and investigating more massive dark-matter dominated 
galaxies may therefore be a more effective way to explore the 
core/cusp issue. These galaxies, fortunately, do exist, and are 
called Low Surface Brightness (LSB) galaxies. 

3.2. LSB Galaxies 

The term LSB galaxies is used here to indicate late-type, 
gas-rich, dark-matter-dominated disk galaxies. Their opti- 
cal component is well-described by an exponential disk with 
an (extrapolated) inclination-corrected central surface bright- 
ness f ainter than no n ~ 23 ma g arcsec" 2 (McGaugh & Bothunl 
fl994t McGaug h etai] [l99l IcteBlok et al.l 119951) . Despite 
their low surface brightness, their integrated luminosity is 
a few magnitudes brighter than that of late-type dwarf 
galaxies (Mb ~ -18 to ~ -20 for LSB galaxies, as op- 
posed to Mb > —16 for the dwar f galaxies). As noted , 
they are gas-rich (Mhi/Lr > 1: ISchombert et alj 119921: 
iMcGaugh & de Blokl[T997t ISchombert et al. 200 1]), and fheir 
interstellar medium has a low metallicity (IMcGaugh 1 1994t 
Ide Blok & van der Hu lst 1998). Their optical appearance is 
dominated by an exponential disk with a young, blue popula- 
tion, with little evidence for a dominant old population. Addi- 
tionally, these galaxies do not have large dominant bulges, and 
seem to have h ad a star formation history with only sporadic 
star formation ([van der Hulst et alJfl993t Ivan den Hoek et all 
I2000t iGerritsen & de Bloklll99'9l Central light concentra- 
tions, if present at all, tend to be only fractionally brighter 
than that of the extrapolated exponential disk. In terms of 
their spatial distribution, they ar e found on the outskirts of the 
large scale structure filaments dBothun et al.ll 1993k iMo et aTl 
1994). In short, most observational evidence indicates that 
these galaxies have had a quiescent evolution, with little evi- 
dence for major merging episodes, interactions, or other pro- 
cesses th at might have stirred up the baryonic and da rk matter 
(see also lBothun et al.lll997tllmpev & Bothunll 19971 

As for the term "LSB galaxies", there is some confu- 
sion in the literature about what type of galaxies it applies 
to. The type of LSB galaxies most commonly studied, 
in particular with regards to the core/cusp controversy, are 
the late-type LSB galaxies whose properties are described 
above. The other type of LSB galaxies often discussed in 
the literature are the massive, early-type, bulge-dominated 



LSB galaxies. These galaxies have p roperties entirely differ- 
ent from the late-typ e LSB galaxies (Sprav berrv et alJ[T995b 
iPickering etaEfT997h . The massive LSB galaxies are a lot 
more luminous and their optical appearance is dominated by 
a brig ht central bulge with a c learly detectable old popula- 
tion (jBeijersbergen et al.l 119991). Man y of them have low- 
level AGN activity (Schombert 1998). All indications are 
that the evolution of these galaxies has been entirely differ- 
ent from that of late-type LSB galaxies: if anything, they re- 
semble SO galaxies with extended disks, rather than late-type 
galaxies. The presence of the dominant bulge also indicates 
that their central dynamics are likely to be dominated by the 
stars, rather than dark matter. In the following, the term "LSB 
galaxies" therefore refers to late-type LSB galaxies only. 

3.3. Early HI observations of LSB galaxies 

The first detailed studies of large samples of LSB galax- 
ies soon led to the picture of them being unevolved, gas- 
rich disk galaxies, as described above. The observation that 
they followed the sam e Tully -Fisher relation as normal galax- 
ies dZwaan et al] 1 1995b was intriguing, as this implied they 
had to be dark-matter dominated. Follow- up radio synthe- 
sis observations in HI (Ide Blok et al.lll996l) soon confirmed 
this. Though the resolution of these early observations was 
limited, the derived rotation curves clearly resembled those 
of late-type dwarf and "normal" disk galaxies: a slow rise, 
followed by a gradual flattening. When expressed in terms 
of scale lengths, the rotation curves of LSB and HSB galax- 
ies of equal luminosity turned out to be very similar, indi- 
cating that LSB galaxies ar e in general low density objects 
(Ide Blok & McGaugh|[T996l) . 

Mass models derived using the rotation curves clearly 
showed that for reasonable assumptions for the stellar mass- 
to-light ratio, T*, the dynam ics of LSB galaxies had to be 
dominated by dark matter (de Blok & McGaugh 1997). As- 
suming that the stars had to dominate the dynamics in the in- 
ner parts (the so-called maximum disk solution) led to unreal- 
istically high T* values, and, even when taken at face value, 
still showed a need for a moderate amount of d ark matter at 
small radii (see also McGaug h & de Bloklll998l) . 

The distribution of the dark matter at fi rst sight seemed sim- 
ilar to that in gas-rich dwarf galaxies (Ide Blok & McG augh 
1 19971). Because of the l imited resolution of the data, 
Ide Blok &~M cGaughi d 19971) did not attempt fits with the 
NFW model, but noted that the halos had to be extended, dif- 
fuse and low density. 

A first attempt at comparing the HI data with CDM predic- 
tions was made by [McGaugh & de Blokj (119981) . Rather than 
making fits to the rotation curve, they simply assumed that the 
typical velocity V200 of the halo had to equal the outer (maxi- 
mum) rotation velocity of the galaxy. The strict cosmological 
relation between c and V200 then automatically yields a value 
of c compatible with ACDM. Adopting these values, the re- 
sulting halo rotation curve turned out to b e very differen t from 
the observed curve, in a similar way as the Moore ( 1994]) anal- 
ysis: the NFW curve is too steep and rises too quickly in the 
inner parts. The only way the halo curve could be made to 
resemble the observed curve, was by abandoning the cosmo- 
logical (c, V200) relation. 

Similar conclusions were derived by ICote et al] (2000). 
They presented high-resolution HI observations of dwarfs in 
the nearby Centaurus and Sculptor groups, and noted that the 
derived rotation curves did not agree with the NFW model. 

A possible explanation that was soon put forward was 
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that there were still unrecognized systematic effects in the 
data, that would give the false impression of a core-like be- 
haviour. Initi a lly, at tention focussed on the resolution of the 
Ide Blok et al.l (119961) HI observations. These had beam sizes 
of ~ 15", resulting in the HI disks of the LSB galaxies in- 
vestigated having a diameter of between 3 and 18 indepen- 
dent beams. This limited resolution can potentially affect the 
shapes of the rotation curves thr ough a proces s called "beam 
smearing", as also mentioned in de Blo k et al.l 0996). In ob- 
servations with limited resolution, the beam smearing process 
decreases the observed velocities (compared to the true veloc- 
ities), and in extreme cases can turn any steeply rising rotation 
curve into a slowly rising solid-body one. This would there- 
fore give the impression of a core being present in the data, 
whi le the true distribution could still be cuspy. In their pa- 
per, Ide Blok & McGaugE (119971) argued, through modelling 
of these beam smearing effects, as well as the direct detec- 
tions of steeply rising rotation curves in the data, that while 
some beam smearing was indeed present, the effect was not 
strong enough to completely "hide" the dynamical signature 
of a cusp, and concluded that the data were consistent with 
the existence of dark matter cores. 

An alternative in terpr etation was given in 
van den Bosch et al.l (|2000l) who used the 
de Blok & M cGaugh d 19971) data, along with high-resolution 
literature rotation curves of a number of late-type "normal" 
and dwarf galaxies, such as DD0154 and NGC 247. They 
derived and applied explicit analytical corrections for beam 
smearing and concluded that the LSB galaxy HI data were 
consistent with both cored and cuspy halos. In their analysis 
of the other gas-rich dwarf and late-type galaxies, they found 
evidence for cores in the dwarf galaxies, but detected a steep 
mass-density slope — consistent with a cusp — in NGC 247, 
the one late-type disk galaxy that met their sample selection 
criteria. 

In a related paper, Ivan den Bosch & Swatersl (120011) derive 
similar conclusions for a different, larger sample of dwarf 
galaxies observed in HI as part of the We sterbork HI Survey 
of Irre gular and Spiral Galaxies (WHISP) (Ivan der Hulst et al] 
120011) . They attempt to correct for adiabatic contraction and 
resolution and conclude that, while in a majority of the galax- 
ies they investigate a core is somewhat preferred in terms of 
fit quality, they cannot exclude halo models that have a steep 
inner mass-density slope. 

Clearly, with this wide range of (sometimes contradictory) 
conclusions, obtaining higher resolution data is the only way 
to put stronger constraints on the exact distribution of dark 
matter in galaxies. 

3.4. Ha observations 

After the initial HI observations, the most efficient way to 
improve the resolution was by obtaining observations in the 
Ha line. These usually took the form of long-slit spectra taken 
along the major axes of the galaxies, and resulted in an order 
of magnitude improvement in resolution. Typical observa- 
tions now had resolutions of a few arcseconds, instead of a 
few tens of arcseconds. 

Some of the first H a long-slit rotat i on cu rves of LSB galax- 
ies were presented by IS waters et al.l (120001) . They found that 
the Ha curves indeed did rise somewhat more steeply in the 
inner parts than the HI curves, and for one or two galaxies 
very much so, but they noted that on the whole the HI and Ha 
curves were consistent when the beam smearing effects were 
taken into account. The Ha curves and corresponding mass 



models did still indicate that LSB galaxies had to be domi- 
nated by dark matter. The curves also still rose less steeply 
than those of higher surface brightness late-type galaxies of a 
similar luminosity, as indicated by the agreement between the 
curves when they were radially scaled using their exponential 
disk scale lengths. 

A large set of high-resolution, Ha long-slit rotatio n 
curves was published an d an alysed bvlMcGaugh et al.l (12001b . 
Ide Blok et all (l2001albl) and ldeBlok & Bosmal d2002l). Many 
of the se galaxies had been part of the de Blok & McGaugh 
( 1997) sample, and a direct comparison between the Ha and 
HI rotation curves showed that in the majority of cases, beam 
smearing effects were present, but not significant enough to 
alter the previous conclusions regarding the dark matter dis- 
tribution. In the few cases where there were large differ- 
ences this could be explained by inclination effects or ellip- 
tical beam shapes in the HI observations. 

The Ha curves thus showed that the inner, slowly ris- 
ing slopes could not be caused by resolution effects. This 
meant that for reasonable (stellar population synthesis in- 
spired) values, LSB galaxies were still dark matter dom- 
inated throughout their disk, even though the Ha curves for- 
mally allowed maximum disk fits with even higher T* values 
than the HI observations did. These high values are, how- 
ever, completely at odds with the observed colours and star 
formation histories, and still imply a significant dark matter 
fraction. 

Comparison of mass models assuming PI and NFW models 
showed that the data were better described by the PI models. 
In many cases, NFW models did yield reasonable fits, but usu- 
ally with very low concentrations and high V200 values, incon- 
sistent with the cosmological (c, V200) relation. Taken at face- 
value these results would imply that the halos of LSB galaxies 
have barely collapsed, and have typical velocities many times 
higher than those of the galaxies that inhabit them. The low 
c-values are, however, due to the intrinsically different shape 
of the NFW rotation curves compared to the solid-body ob- 
served curves. An NFW curve has an inner velocity slope 
v oc r 1 / 2 , and the only way this kind of curve can fit the ob- 
served solid-body v oc r curve is by stretching, resulting in the 
low c and high V200 values. 

The observational inne r mass density slopes were derived 
by Ide Blok et al. (2001a) and plotted against the resolution 
(physical radius of the innermost point) of the rotation curve. 
A comparison with the expected slopes from various halo 
models showed that the majority of the data scattered towards 
the predicted slopes of the PI model. They also showed that 
for resolutions of ~ 1 kpc the PI and NFW models yielded 
identical slopes. This would go some way towards explain- 
ing why some of the lower resolution HI observations were 
unable to distinguish between models. The behaviour of 
the change in slope when going from the lower resolution 
HI observation to higher resolution Ha observations is con- 
sistent with that expecte d from the PI model, as shown in 
Ide Blo k & Bosma (2002). Clearly, for unambiguous determi- 
nations of the inner mass density slope, resolutions of better 
than a kpc are needed. 

Independ ent observations and an alyses came to similar con- 
clusions. Marchesini et all (|2002|) obtained longslit H a ob- 
servations of some of the Ide Blok & McGaugh ! (119971) galax- 
ies, and also found little evidence for strong beam smearing, 
as we ll as strong evidence fo r the existance of dark matter 
cores. IZackrisson et al.l d2006) obtained long-slit observations 
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of the rotation curves of six extremely blue and bulge-less 
LSB galaxies and also finds strong evidence for the pres- 
enc e of cores in th e dark matter distribution of these galax- 
ies. ISalucci (2001) analysed the rotation curves of over a hun- 
dred rotation curves of disk galaxies, and found clear signa- 
tures for the existence of dark matter cores in these galaxies. 
This analysis is complementary to the work on LSB galax- 
ies, as it also analyses le ss dark-matter-dominated galaxies. 
iBorriello & Salu cci (2001) analysed Ha rotation curves of a 
number of dark matter dominated galaxies, and also found 
strong evidence for core-like dark matter distributions. They 
also fit Burkert haloes to their data and find that these also 
provide good fits. Analysis for a sample of ~ 25 spiral galax- 
ies is presented in Dona to et alj (|2004), and leads to similar 
conclusions. They also find that the core radius of the dark 
m atter distributi o n is re lated to the disk scale length. 

lHayashi et al.l (120041) find th at the Ha rotat ion curves are 
consistent with steep slopes. Ide Blokl (120051) suggests that 
these conclusions are based on artificial constraints imposed 
on the fitting functions. Once these constraints are removed, 
17 of the 20 galaxies with the highest-quality data are best 
fitted with cored models. Three are possibly consistent with 
cuspy models; two of these are high-surface brightness dwarf 
galaxies that are likely dominated by stars. 

3.5. Possible systematic effects 

At first glance, the observational data seems to provide 
good evidence for the presence of cores in LSB galaxies. This 
does of course not exclude the possibility that systematic ef- 
fects might be present in the data that could give the (false) 
impression of cores. A number of studies have therefore fo- 
cused on these effects and asked whether the detected cores 
could still be consistent with cuspy halo models. 

The systematic effects that were investigated fall into two 
categories. 

1 . Pointing problems — If small offsets exist between the 
central slit position and the true, dynamical center, this 
will cause the spectrograph slit to miss the cusp. This 
could be due to inaccurate telescope pointings or, alter- 
natively, if large physical offsets between the dynamical 
and photometric centers of galaxies exist, the measured 
slopes will also be biased downwards. 

2. Non-circular motions — The fundamental assumption 
in the observational analyses discussed so far is that the 
gas moves on circular orbits. If for some reason the 
orbits are elliptical, or of the gas motions are disturbed, 
this will also lead to an underestimate of the slope. 

These effects are difficult to recognize in long-slit, one- 
dimensional Ha data without additional information. Many 
authors therefore, apart from modeling these effects, also em- 
phasized the need for high-resolution, two-dimensional ve- 
locity fields. These make pointing problems irrelevant, while 
non-circular motions can be directly measured. Fortunately, 
these velocity fields are now available (see Sect. 3.6), largely 
superseding the results derived from the Ha rotation curves. 
Nevertheless, for completeness, the further analysis of the Ha 
curv es is briefly di s cusse d here. 

In Ide Blok et al.l d2003) a first attempt was made at model- 
ing the observational systematic effects. Their main conclu- 
sion was that NFW halos can be made to resemble dark matter 
cores only if, either systematic non-circular motions with an 



amplitude of ~ 20 km s" 1 exist in all disks, or systematic tele- 
scope pointing offsets of ^3-4" exist for all observations, or 
if the dynamical and photometric centers are systematically 

o ffset in all galaxie s by ~ .5 - 1 kpc. 

iMarchesini et al.l (120021) and Ide Blok & Bosmal (120021) 
compare independent sets of observations of the same galax- 
ies, obtained at different telescopes, by independent groups, 
using independent data sets, and find no evidence for tele- 
scope pointing errors. In general, galaxies can be acquired 
and positioned on the slit with a repeatability accuracy of 0.3" 
or so. 

Using further modeling, deBloketal. (2003) find that a 
halo model with a mildly cuspy slope a = -0.2 ± 0.2 gives the 
best description of the data in the presence of realistic obser- 
vational effects. It is interesting to note that this best-fit ting 
slope had alre ady been derived bvlKravtsov etaf] (119981) on 
the basis of the Ide Blok et al J (119961) dat a. 

An analysis by Spekke ns et al.l ([2005) of longslit Ha rota- 
tion curves of 165 low-mass galaxies comes to the same con- 
clusion: depending on how they select their sample, they find 
best-fitting slopes of a = -0.22 ± 0.08 to a = -0.28 ± 0.06. 
They also model pointing and slit offsets, but come to the 
conclusion that, after correcti o n, the ir data are consistent with 
cuspy haloes. Swaters et al. (2003a) also present extensive 
modeling of high-resolution long-slit Ha rotation curves, and 
show that while their data are consistent with a = cores, 
steeper slopes cannot be ruled out. 

Given the difference in interpretation of otherwise similar 
samples, a double-blind analysis of the modeling performed 
by the various groups would have been interesting. However, 
with the availability of high-resolution velocity fields, this is 
n ow a moot i s sue. 

Rhee et al] (120041) attempt to model many of the observa- 
tional effects using numerical simulations. Most of their con- 
clusions apply to long-slit observations. The systematic ef- 
fects they investigate (dealing with inclination effects, non- 
circular motions and profile shapes) have now been directly 
tested on high-resolutio n velocity fields and do not criti- 
cally affect the data (e.g. .[Gentile et al.l2 005; Zackrisso n et al. 
1 20061: iTrachternach et alJl2008HOh et alJl2008t Ide Blok et al. 
l2008tlKuzio de Narav et alJl2008Ll2009l) . 

3.6. High- resolution velocity fields 

As noted, high-resolution two-dimensional velocity fields 
provide the context which long-slit observations are missing. 
With these velocity fields the pointing problem becomes irrel- 
evant, offsets between kinematical and dynamical centers can 
be directly measured, as can non-circular motions. Following 
is a brief overview of the various observational studies that 
have presented and analysed these velocity fields within the 
context of the core/cusp debate. 

Some of the first high-resolution optical velocity fields 
of late-type dwarf a nd L SB galaxies were presented by 
iBlais-Ouellette et all (120011) . They analyse Ha Fabry-Perot 
data of IC 2574 and NGC 3109, two nearby dwarf galax- 
ies, and derive slowly rising rotation cu rves, consistent with 
a core . This work was later expanded in Bla is-Ouellette et al] 
(2004), and led to the work presented in Spano et al. (2008), 
where optical velocity fields of 36 galaxies of different mor- 
phological types are presented. All three studies find that the 
PI model generally provided better fits than NFW models. If 
NFW models give fits of comparable quality, then this is usu- 
ally at the cost of an unrealistically low value and non- 
cosmological (c,V2oo) values. A similarly large collection of 
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Ha velocity fields is presented in Dica ire et al.l (J2008). They 
do not explicitly address the core/cusp issue, but show that 
when bars are present, their influence on the velocity fields is 
v ery noticable. 

iKuzio de Narav et alj (120081 120091) present DensePak ve- 
locity fields of LSB galaxie s, many of them taken from the 
Ide Blok &~ McGaugh d 19971) sample. Their conclusions are 
that NFW models provide a worse fit than PI models, for all 
values of Y*. Where an NFW model could be fit, the c-values 
generally again do not match the cosmological CDM (c, V200) 
relation. They introduce a "cusp mass excess": when the pre- 
dicted (c, V200) relation is assumed, and the V200 velocities are 
matched with those in the outermost observed velocities, the 
inner parts require about ~ 2 times more dark matter mass 
th an is implied by the ob s erved r otation curves. 

IKuzio de Naray et all (120081 120091) also explore non- 
circular motions: they find that random velocities with an 
amplitude of ~ 20 km s" 1 are needed to bring the observed 
curves in agreement with the CDM predictions. A comparison 
of simulated long-slit observations (extracted from the veloc- 
ity fiel ds) with the original long-slit data from McGaug h et al.l 
d2001l) and Ide Blok & Bosnia! (120021) shows good agreement. 

Swaters et al. (2003b) present a DensePak velocity field of 
the late-type dwarf galaxy DDO 39. They derive a rotation 
curve, and show that its slope is steeper than implied by lower 
resolution HI data, and also differ ent from earlier long-slit 
data from lde Blok & Bosmal (|2002). They indicate that mea- 
surable non-circular motions are present, but do not explicitly 
quantify them. They show fit results for NFW halos, but do 
not show the corresponding results for a core model, so fur- 
th er comparisons are dif ficult to make. 

IWeldrake etaTJ (120031) present the HI velocity field of NGC 
6822, a nearby Local Group dwarf galaxy. Their data has a 
linear resolution of about 20 pc, so beam smearing is defi- 
nitely not a problem. The rotation curve shows a strong pref- 
erence for a core-like mo del, but they do not quantify possi- 
ble non-circular motions. ISalucci et al.l (120031) present simi- 
larly high-resolution HI data (including VLA B array) of the 
dwarf galaxy DDO 47. Their analysis shows the dynamics 
are not consi s tent w ith a cusp. Similar results are derived in 
iGentile et al.l (12 004) for a number of spiral galaxies. 

Sim on et al.1 (|2005) present results from CO and Ha veloc- 
ity f ields of 5 low-mas s da rk-matter-dominated galaxies (see 
also ISimon e"taill2003l and iBolatto et al.ll2002l) . They derive 
a range of slopes, from core-like (a = 0.01 for NGC 2976) 
to cuspy (a = 1.20 for NGC 5963). Note that NGC 5963 has 
an inner bright disk, and it may therefore not be dark-matter- 
dominated all the way to the center, ma king the value of it s 
mass-density slope uncertain (see also iBosma et al.l 119881) . 
The average slope they derive is a = 0.73 ±0.44. Their analy- 
sis method differs in a few aspects from the other studies ref- 
erenced in this paper. Firstly, their inner slope values are de- 
rived from a single power-law fit to the entire rotation curve. 
Most of their models do not take into account the gas compo- 
nent due to a lack of HI observations. In low-mass galaxies 
the gas can dynamically be more important than the stars (es- 
pecially in the outer parts), and correcting for this component 
could potentially change the derived slopes. 

More importantly, the mass models in Simon et al. (2005) 
are derived under the explicit assumption of a constant in- 
clination and position angle for each galaxy. Most velocity 
fields of nearby galaxies show radial inclinatio n and position 
angle trends, especially in the outer parts (e.g.. Ide Blok et all 



120081) . ISimon et all (120051) attribute these to radial velocities 
that give the impression of changes in inclination and position 
angle. Nevertheless, the non-circular velocities derived in this 
way are typically less than ~ 20 km s" 1 , and in a few galax- 
ies less than ~ 5 km s" 1 . Whether these velocities are real or 
whether inclination and position angle changes are preferred 
would be an interesting topic of further study. Simo n et al.1 
(2005) also derive harmonic decompositions of the velocity 
fields to study non-circular motions, but they do not list the 
values of the harmonic coefficients. It would be similarly in- 
teresting to compare these with results derived for other disk 
and LSB galaxies. 

Direct measuremen ts of non-circular motions are presented 
in IGentile et al.l (120051) . Using high-resolution HI observa- 
tions, they make a harmonic decomposition of the velocity 
field, and show that non-circular motions are only present at 
the level of a few km s" 1 . This is a factor of ~ 10 lower than 
is needed for the non-circular motions to wipe out the kine- 
matical signature of a cusp and to give the impression of a 
core. 

An analysis by IGentile et al.l ((2007a) of another gas-rich 
dwarf galaxy, NGC 3741, based on the entire 3D HI data 
cube, showed non-circular motions around 5-10 km s" 1 , and 
a strong preference for a core model. NFW models could be 
accommodated, but with the usual caveat of the fit parameters 
n ot being consistent with the cosmological (c, V200) relation. 

iTrac hternach et al.1 (120081) present harmonic decomposi- 
tions of the velocity fields of galaxies from T he HI 
Nearby Galaxy Survey (THINGS: Walter et al.l 12008). The 
THINGS survey covers a large range in galaxy properties, 
from luminous ear l y-type disk galaxies to late-type dwarfs. 
ITrachternach et all (|2()()8 ) find a relation between the median 
strength of the non-circular motions and the luminosity of the 
galaxies, indicating that the non-circular motions are asso- 
ciated with the baryons. Luminous disk galaxies have non- 
circular motions up to ~ 30 km s , mostly associated with 
bars and spiral arms. These then decrease rapidly to a level of 
only a few km s" 1 for dwarf galaxies like DDO 154. They also 
found that offsets between photometric and kinematic centers 
were typically ~ 200 pc or less. 

This low level of n on-circular motions seems inconsistent 
with observations by iPizzella et al.l (120081) . who find signif- 
icant non-circular motions in a sample of 4 LSB galaxies. 
However, their galaxies all contain bright bulges and are 
therefore probably more representative of the class of giant 
LSB galaxies, rather than th e late-type ones disc ussed in this 
paper (in t his regard see al s olCoccato et al.ll2008h. 

Finally, McGaugh etaD (120071) and IGentile et al.l (l2007bl) 
show that the core/cusp issue is not limited to the very in- 
ner parts of galaxies. As mentioned before, the shape of 
the NFW curve is fundamentally different from that of ob- 
served cur ves. One can try and w ork out the implications in 
two ways: iMcGaugh et al.1 (120071) match the observed outer 
rotation velocities with those of corresponding cosmological 
NFW halos (effectively identifying NFW halos that have a 
similar dark matter density at these outer radii as real galax- 
ies), and find that these halos are too massive. Due to the cusp 
mass excess their average density is then also a factor of 2-3 
too high (see also Ide Blok et al.ll2008l) . IGentile et al.1 (12007b!) 
take a different approach and identify halos that have the same 
enclosed mass as observed galaxies. Again, due to the cusp 
mass excess, this results in halos that, compared to what can 
be derived for real galaxies, are more dense in the center and 
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less dense in the outer parts. Due to the fundamentally dif- 
ferent shapes of the mass distributions, the core/cusp issue is 
therefore not limited to the inner parts, but is relevant at all 
radii. 

4. EFFECTS OF BARYONS AND TRIAXIALITY 

The high-resolution velocity fields thus seem to indicate 
mass distributions that are cored. As many observational ef- 
fects such as pointing, center offsets, and even non-circular 
motions seem to be too small to be the cause of the observed 
cores, many studies have attempted to explain the apparent 
presence of cores as the result of processes such as interac- 
tions between and merging of dark matter halos, or the effect 
of baryons on the dark matter distribution. Following is a brief 
discussion of some of these effects. 

4.1. Feedback and merging 

It was already noted that violent, large-scale star formation 
might be able to explain the cores in gas-rich dwarf galaxies, 
but not in more massive galaxies (Navarro et al 1996). The 
implied large bursts of star formation are also not consistent 
with the quiescent evolution of LSB galaxies. An alterna- 
tive way to achieve the r emoval of the cusp was proposed by 
Weinberg & Katzl J2002). They model a rotating rigid bar in 
a disk which is embedded in a CDM halo and show that this 
bar creates a wake in the dark matter. This trailing wake slows 
down the bar, transferring some of the bar angular momentum 
to the dark matter. This "puffs up" the dark matter distribu- 
tio n thus forming a core . Recent numerical simulation work 
by Dubins ki et al.l (|2009) does, however, suggest that the bar 
does not destroy the central cusp, and may even increase the 
halo density slightly. 

A very detailed ca se-study of the effects o f bars and feed- 
back is presented in IValenzuela et alj (120071) . They attempt 
to reproduce the observed rotation curves of NGC 6822 and 
NGC 3 109 using NFW halos and a variety of feedback effects 
and non-circular motions. They present numerical models in- 
cluding gas dynamics, and tune them to resemble the two tar- 
get galaxies. Whilst they succeed in matching the rotation 
curves, it is not clear whether these results can be general- 
ized, due to, e.g., the different assumptions that are made in 
determining inclinations and position angles of the models, 
compared to the data. Although the authors present simulated 
velocity fields, they do not make a direct comparison between 
observed and simulated velocity fields, which would provide 
m uch additional info rmation on the bar and feedback effects. 

Deke l et all (l2003alfbl) argue that merging of cuspy halos in- 
evitably leads to a cusp. The only way to prevent this from 
happening is to puff up the dark matter distribution of the in- 
falling halos before they merge. This way the halos get dis- 
rupted more easily and a core-like distribution can be gradu- 
ally built up. The authors describe the scenario as speculative, 
however, and note that it is unlikely that supernovae will be 
able to cause the puffing up for any galaxy with V > 100 km 
s" 1 . iBovlan-Kolchin & Mai d2004) show that for this process 
to work, none of the merging halos can be cuspy to start with, 
as only mergers between cored halos give a cored merger 
product. Any merg er involving a c usp inevitably leads to a 
cuspy end-product. Dehnen (2005) also shows that the final 
slope always equals that of the steepest component. Setting up 
cored halos and ensuring that they remain cored is apparently 
not trivial. 

4.2. Dynamical Friction 



El-Zantetd]([2001) propose a different way to make ha- 
los cored: they note that merging gas clouds of ~ 10 5 M Q 
(for dwarfs) to ~ 10 8 M Q (for spirals) can disrupt cusps 
through dynamical friction. If this happens early enough in 
the universe (when halos were smaller), this proc ess could be 
very e ffi cient. Similar scena r ios are presented in Toni ni et alj 
(2006). iRomano-Dfaz et alj d2008l) in their study also argue 
that the effect of the baryons on the halo structure must be 
significant. They suggest that in the presence of baryons, ini- 
tially a very steep cusp is formed (with a ~ -2), which is then 
heated by sub-halos through dynamical friction, and, subse- 
quently, from the inside out becomes shallower. According 
to their analysis, the end result is a density profile that is less 
steep than a = -1 in the inner few kpc, and may even be cored 
in the very center. 

iJardel & Sellwoo J (120091) take the opposite view. They 
model the dynamical friction process, but argue that it is dif- 
ficult to find "free-floating" baryon clumps massive enough 
to make this process happen, as these clumps must not be as- 
sociated with dark matter (due to the cusp that will then be 
forme d; see above). They note that work by Kaufmann et al. 
(2006) implies clump masses that are two orders of magnitude 
too small. 

The proces ses just describe d are simila r to those 
proposed by iMashchenko et aT] d2006l [2008); see also 
iRicotti & Wilkinson! d2004l) . Thev make it all happen in the 
early universe, when mass scales and the required amount 
of baryons were both smaller. Their numerical simulations 
suggest that cusps can be erased in the very early universe 
(z < 10) when (proto-) galaxies had approximately the size of 
the HI holes and shells observed in the disks of present-day, 
gas-rich spiral and dwarf galaxies. In objects that small, a 
random motion of only ~ 7 km s" 1 (i.e., equal to the HI ve- 
locity dispersion observed in local disk galaxies) is sufficient 
to disrupt the cusp and keep the halos co red until the present 
day. N ote though that in a recent analysis. ICeverino & K lypin 
( 2009) perform similar calculations, but find that halos remain 
cuspy, and suggest that differences in the simulation approach 
m ight be responsible for thi s. 

Chen & McGaugh (2008) also argue against the idea of cre- 
ating cores at high redshift. Their argument is that cuspy halos 
cannot explain gravitational lensing results, but demand halos 
with an even steeper mass distribution (a singular isothermal 
sphere with slope a = -2). This means that if the halos of the 
elliptical galaxies that do the lensing indeed form with NFW- 
like profiles, these need to steepen over the course of their 
evolution, using some process for which the quiescent adia- 
batic contraction is the most likely candidate. LSB galaxies, 
on the other hand, need to experience vigorous bursts of early 
star formation that drive feedback to erase the initial cusp. 
These scenarios are at odds with what we know about the evo- 
lution of these galaxies: ellipticals typically undergo a large 
amount of merging, with the associated vigorous star forma- 
tion, as evidenced by their extensive old stellar populations. 
LSB galaxies show no evidence at all for any kind of violent 
interactions, nor intense star formation. Dynamical friction 
and adiabatic contraction thus seem to place demands on the 
evolution of elliptical and LSB galaxies contrary to what can 
be derived from their star formation histories. 

4.3. Triaxiality 

lHavashi & Navarro! d2006l) and lHavashi et alj d2007l) show 
that introducing an elliptical disturbance in an NFW poten- 
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tial can lead to systematic (non-circular) motions that, when 
unrecognized, could be interpreted as evidence for a core. 
Their arguments particularly apply to longslit Ha observa- 
tions, where indeed the context of a velocity field is not avail- 
able to gauge the validity of the circular motion assump- 
tion. As described earlier, work on high-resolution vel ocity 
fields bvlGentile et all (12005b . iTrachternach et al.l (120081) and 
iKuzio de Narav et al.l d2009t) has put significant observational 
constraints on the strength of non-circular motions and the el- 
lipticity of the (equatorial) potential. For a large sample of 
disk and dwarf galaxies, this potential is consistent with be- 
ing round and the non-circular motions are too small to give 
the illusion of cores in lo ng-slit observations. 

The resul ts fro m lHavashi & Navarrol d2006l) and 
lHavashi et alj (2007) assume massless disks. iBailin et all 
(120071) present results from simulations of self-consistent 
massive disks in triaxial halos and find that the baryons 
circularize the potential rapidly, even for low-mas s disks, 
thus w iping out a large part of the tri-axiality signal. IWidrowl 
(2008) also presents models of "live" disks in tri-axial halos, 
and can approximate an observed LSB rotation curve by 
introducing tri-axiality in the halo. 

In all these studies it would be interesting to compare sim- 
ulated velocity fields with the observed ones. As has become 
clear from the preceding discussion, modeling the long-slit 
rotation curves leaves too many ambiguities which only stud- 
ies of the velocity fields can addr ess. In an observational 
study, Kuz io de Narav et all Q2009) subject model NFW ve- 
locity fields to the DensePak observing procedure, and show 
that even in the presence of observational uncertainties the 
signature of an NFW velocity field can be observed. 

They show that axisymmetric NFW velocity fields are un- 
able to reproduce the observed velocity fields for any com- 
bination of inclination and viewing angle. They derive NFW 
velocity fields in an elliptical potential, and show that the only 
way these can be made to resemble the observations is by hav- 
ing the observer's line of sight along the minor axis of the 
potential for 6 out of the 7 galaxies investigated, inconsistent 
w ith a random distributi o n of th e line of sight. 

Kuzio de Narav et al. (2009) also note that the kind of 
rapidl y varying ellip ticity of the potential as proposed by 
lHavashi et alj 0007) might help in better describing the data, 
but that the problem of a preferred viewing angle will remain. 
For an elliptical potential with a random viewing angle, one 
also expects, once in a while, to observe a rotation curve that 
is steeper than the corresponding axisymmetric NFW profile. 



Rotation curves like that are, however, exceedingly rare, if not 
absent, in the available observations. 

5. SUMMARY 

Rotation curves of LSB and late-type, gas-rich dwarf galax- 
ies indicate the presence of constant-density or mildy cuspy 
(a ~ -0.2) dark matter cores, contradicting the predictions of 
cosmological simulations. The most recent simulations still 
indicate resolved mass density slopes that are too steep to be 
easily reconciled with the observations (typically a ~ -0.8 at 
a radius ^0.1 kpc). Claims of shallow slopes at even smaller 
radii depend on the validity of the analytical description cho- 
sen for the mass-density profile. 

Whereas early HI observations and long-slit Ha rotation 
curves still left some room for observational (pointing, res- 
olution) or physical (non-circular motions, tri-axiality) sys- 
tematic effects to create the illusion of cores in the presence 
of a cuspy mass distribution, the high-resolution optical and 
HI velocity fields that have since become available signifi- 
cantly reduce the potential impact of these effects. Measured 
non-circular motions and potential ellipticities are too small 
to create the illusion of a core in an intrinsically cuspy halo. 

This indicates either that halos did not have cusps to be- 
gin with, or that an as yet not understood subtle interplay be- 
tween dark matter and baryons wipes out the cusp, where the 
quiescent evolution of LSB galaxies severely limits the form 
this interplay can take. Adiabatic contraction and dynamical 
friction yield contradictory results, while models of massless 
disks in tri-axial halos result in preferred viewing directions. 
LSB galaxy disks, despite their low T* values, are not en- 
tirely massless, and observations and simulations will need to 
take this into account. Similarly, the difficulties in reconcil- 
ing a possible underlying triaxial potential with the circular- 
izing effects of the baryons also needs to be investigated. In 
short, studies which, constrained and informed by the high- 
quality observations now available, self-consistently describe 
and model the interactions between the dark matter and the 
baryons in a cosmological context are likely the way forward 
in resolving the core/cusp problem. 

I thank the anonymous referees for the constructive com- 
ments. The work of WJGdB is based upon research supported 
by the South African Research Chairs Initiative of the De- 
partment of Science and Technology and National Research 
Foundation. 
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